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ABSTRACT 

Using the idea of smooth crossover from the hadronic matter with hyperons to quark matter with 
strangeness, we show that the maximum mass (M max ) of neutron stars with quark matter core can 
be larger than those without quark matter core. This is in contrast to the conventional softening of 
equation of state due to exotic components at high density. Essential conditions to reach our conclusion 
are (i) the crossover takes place at relatively low densities, around 3 times the normal nuclear density, 
and (ii) the quark matter is strongly interacting in the crossover region. By these, the pressure of the 
system can be greater than that of purely hadronic matter at a given baryon density in the crossover 
density region and leads to M max greater than 2 solar mass. This conclusion is insensitive to the 
different choice of the hadronic equation of state with hyperons. Several implications of this result to 
the nuclear incompressibility, the hyperon mixing, and the neutrino cooling are also remarked. 
Subject headings: dense matter, elementary particles, equation of state, stars: neutron 



1. INTRODUCTION 

Since neutron stars (NSs) are gravitationally bound, 
observation of the NS mass (M) provides a precious 
probe for the equation of state (EOS) and the relevant 
composition and structure of dense matter. The larger M 
requires a stiffer EOS and thereby gives a stringent con- 
dition on the state of matter, since various new phases 
proposed so far, such as pion condensation, kaon con- 
densation, hyperon (Y) mixing, deconfmed quark matter 
and so on, are believed to soften the EOS. 

For example, the EOS is softened dramatically when 
Y takes part in NS cores and the theoretical maximum 
mass fails to exceed the observed mass, M Q b s = 1.44M© 
for PSR1913+16 with M Q being the solar mass. This 
contradiction between theory and observation may in- 
dicate a missing repulsion in hyperonic matter such 
as the 3-body interaction (TNI ) acting universally i n 
the nucleon and hypero n sectors (|Nishizaki et al.l (|2001ft ; 
iTakatsuka et al.l (|2006ll2008H ). 

Very recently, the observation of heavy NS 
with M = (1 .9 7 ± 0.04)M^ for PSR J1614-2230 
(jDemorest et al.l (|2010l )) has been reported, which 
causes even stronger constraint on the EOS. In par- 
ticular, it stimulates intensive discussions on whether 
such a heavy NS can have exotic components in the 
central co r e (lOzel et al.l (120101) : iKurkela et al.l (12010ft: 



iKim et al.l (12011ft: iKlahn et all (12011ft : lLattimer et al. 

(l20rTh: 



(2010); "iWeissenborn et ail 
(12011: IChen et all (1 2 Oil 
iWhittenburv et al.l (j20ll V 



iSchramm et al.l 



Bonanno et al 
^20ll 



In particular, in many of the previous works treat- 
ing the transition from the hadron phase to the quark 



phase at high density, point-like hadrons and quarks are 
treated as independent degrees of freedom, and Gibbs 
phase equilibrium conditions are imposed. However, 
quantum chromodynamics (QCD) tells us that such a 
description is not fully justified, since all hadrons are ex- 
tended objects composed of quarks and gluons. Indeed, 
one may expect a gradual onset of quark degrees of free- 
dom in dense matter associated with the percolation of 
finite-size hadrons, i.e., a smooth crossover from hadronic 
matter to q uark ma t ter: S e e seminal works o n hadron 
percolation ([Bavm I (| 1979ft : iCelik et al.l (11980ft V recen t 
works on hadron-quark crossover (jBavm et all ((2008ft : 
iMaeda et al.l ([2009)) an d a recent review on p hase tran- 
sitions in dense QCD (Fukushi ma et al.l (|2011ft ). 

The aim of this paper is to make a phenomenological 
study on the mechanism realizing hybrid stars (NS with 
a quark core) compatible with 2Mq from the point of 
view of the smooth crossover between the hadronic mat- 
ter with hyperons and the quark matter with strangeness. 
Unlike the traditional approach to the hadron-quark 
transition, the crossover picture gives us a novel tool 
to study whether the hybrid star with a quark core is 
compa tible with massive neutron stars ([Takatsuka et al.l 

2. HADRONIC MATTER WITH STRANGENESS 

Let us start with an example of hadronic EO S 
with hyperons (H-EO S) (|Nishizaki et al.l ([20011 l200l : 
ITakatsuka et al.l ([2006ft ) obtained with the following pro- 
cedure: (i) Effective two-baryon potentials Vbb> (B = n, 
p, A, £~) are constructed on the basis of the G-matrix 
formalism to take into account their density-dependence, 
(ii) A phenomenological 3-body nucleon force Unn 1 ex- 
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presse d in a form of two-body potential (jFriedman et al.l 
(|198lD ) is introduced to reproduce the saturation of sym- 
metric nuclear matter (the saturation density po = 0.17 
fm~ 3 and the binding energy Eq = —16 MeV) and the in- 
compressibility n. (iii) We assume universal three-body 
repulsion even for the hyperons through the replacement, 

Unn' -^Ubb 1 - (iv) By using Vbb 1 + Ubb', we obtain the 
hadronic EOS under charge neutrality and /3-equilibrium 
and calculate particle composition j/j (i — n, p, A, E — , 
e~ and pT ) as a function of total baryon density p. This 
gives our H-EOS of neutron star matter with hyperon 
mixing. 

The stiffness of H-EOS is specified by k; the case with 
K = 300 (250) MeV is denoted by TNI3u (TNI2u). 
Here the subscript "u" means that the 3-body force 
is introduced universally including the hyperon sector. 
The maximum mass (M max ) of the NSs with hyperon- 
mixed core with universal 3-body force becomes M max ~ 
1.82(1.52)M Q for TNI3u (TNI2u): On the other hand, 
with only the nucleon 3-body force (TNI3), it reduces 
to Mmax — I.IMq which is not even compatible with 
Mobs = 1.44M . One of the effects of the universal 
3-body force is to delay the onset of hyperon mixing, 
e.g. pa ~ p-%- ~ Ap for TNI3u and TNI2u cases, while 
Ps- =s 2.2p , p A ~ 2.5p for TNI3. 

To check the mo del- dependence of the H-EOS, we con- 
sider two alternat ive equation s of st ate with hyperons: 
AV18+TBF+AE (jBaldo et al.l ((2000)) which is based on 
the G-matrix approach with the AV18 nucleon-nucleon 
potential, Urbana-type three-body nucleon potential 
and the Nijmegen soft-core nucleo n- hyperon potential; 
SCL3AE (|Tsubakihara et al.l pOlOh ) which is a relativis- 
tic mean-field model with chiral SU(3) symmetry. The 
maximum masses (central baryon densities, radius) of 
NSs composed only by hadronic EOSs, AV18+TBF+AE 
and SCL3AE, are M max - 1.22M Q (7.35p , 10.46km) 
and M max ~ 1.36M (5.89p , 11.42km), respectively. 

In Fig.l, we plot the energy per particle (E/A) of 
hadronic EOS with hyperons (TNI3u, TNI2u, TNI3, and 
SCL3AE) toge ther with the nucle ar EOS without hy- 
perons, APR (|Akmal et al.l (|1998f )). We do not show 
AV18+TBF+AE since it is almost the same with TNI3. 
Filled circles on each line denote the density where the 
hyperons start to mix. From this figure, one can see that 
(i) the mixture of hyperons softens the equation of state 
relative to APR, and (ii) onset of the hyperon mixture 
is shifted to higher density if we consider the universal 
3-body interaction. In the following, we adopt TNI3u 
as a typical example of the H-EOS and comment on the 
results with other H-EOSs at the end. 

3. QUARK MATTER WITH STRANGENESS 

Next we consider EOS of quark matter with 
strangeness (Q-EOS). Around the baryon density of a 
few times po where the hadron-quark crossover is ex- 
pected to take place at zero temperature, the deconfincd 
(or partially deconfined) quarks would be still strongly 
interacting. Analogous situation at high temperature has 
been expected theoretically and evidences of a strongly 
interacting quark-gluon plasma (sQGP) were found ex- 
periment ally in relativistic heavy -ion collisions (see e.g. 
a review, iFukushima eTail (|2l)Tlt )L 

Since numerical simulations of quantum chromody- 
namics (QCD) on the space-time lattice at high baryon 




p/p 





Fig. 1. — Energy per particle (E/A) as a function of the to- 
tal baryon density p. Solid lines denote EOSs with hyperon mix- 
ing; TNI3u (universal 3-baryon force with re = 300MeV), TNI3 (3- 
nucleon force with re = 300MeV), TNI2u (universal 3-baryon force 
with re = 250MeV), SCL3AS (relativistic mean-field model with 
chiral SU(3) symmetry). Dotted lines denote the EOS without hy- 
peron, APR (AV18+5u+ UIX*). Filled circles on each line show 
the density where the hyperons start to mix. 

density is not available due to the notorious sign 
problem, we treat the strongly interacting quark mat- 
ter (sQM) at zero temperature by the (2+l)-flavor 
Nambu— Jon a -Lasin i o (NJL) m o del (s e e the reviews, 
(iVogel et all (1199111; iKlevanskvl (j 19921) ; lHatsuda et all 
(1994); Buballaf|2005j)). It is an effective theory of QCD 
and is particularly useful for taking into account the 
non-perturbative phenomena such as the partial restora- 
tion of chiral symmetry at high density. The model La- 
grangian reads 

£n.tl =q(i0 - m)q + ^ ^[(?A>) 2 + (qij 5 \ a q) 2 } 

+G D [dctq(l + l5 )q + h.c] - ^-(qrq) 2 , (1) 

where the quark field qi (i — u, d, s) has three colors 
and three flavors with the current quark mass m^. The 
term proportional to Gs is a !7(3)l x U(3)r symmetric 
four-fermi interacti on w here A a are the Gell-Mann ma- 
trices with A = y/2/3I. The term proportional to Gc 
is the Kobayashi— Maskawa— 't Hooft (KMT) six-fermi 
interaction which breaks £^(1)^ symmetry. The third 
term proportional to g v is a phenomenological vector- 
type interaction. It has some varieties depending on its 
flavor-structure: Here we use the form given in Eq.([T]) 
which leads to an universal flavor-independent repulsion 
among quarks. 

In the mean-field approximation, the constituent quark 
masses Mi (i = u, rf, s) arc generated dynamically 
through the NJL interactions (Gs,d), Mi = mi — 2Gs<Ji — 
2GD&i<Jk, where Oi — (q~iqi) is the quark condensate in 
each flavor, and (z, j, k) corresponds to the cyclic per- 
mutation of u, d and s. On the other hand, the vector 
interaction (g v ) leads to an effective chemical potential 

(|Asakawa et al.l (fl989Tn . fi\ eS) =^-g v £i («<«>• 

Basic parameters of the NJL model are determined 
from hadron phenomenology : In this paper, we ad opt so- 
called the HK parameter set ([Hatsuda et alj (|T994)), A = 
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631.4 MeV, G S A 2 = 1.835, G D A 5 = 9.29, m u . d = 5.5 
MeV, m s = 135.7 MeV, where A is the three-momentum 
cutoff. The magnitude of g v has not been well deter- 
mined, but recent studies of the model applied to the 
QCD phase diagram suggest th at it can be compara - 
ble to or even l arger than Gs (jBratovic et al.l (2012); 
lLourenco et al.l (|2012fl ). so that we change its value in 
the range, 

0<|^<1.5. (2) 
Gs 

The EOS of strongly interacting quark matter with 
strangeness is obtained from the above model under 
charge neutrality and /3-equilibrium with u, d, s, e~ and 
fi~ . It turns out that the strange quarks starts to appear 
at p ~ 4/9o regardless of the magnitude of g v . Also, pT 
never appears under the presence of the strangeness. 

4. HADRON-QUARK CROSSOVER 

To realize the smooth transition between the hadronic 
matter and the quark matter, we first define a "crossover 
density region" characterized by its central value p and 
its width r. The description of the matter in terms of 
the pure hadronic EOS is accurate for pCp-T, while 
the description in terms of pure quark EOS is accurate 

for p ^> p + r. In the crossover region, p — T ^ p^ p + 
r, both hadrons and quarks are strongly interacting, so 
that neither pure hadronic EOS nor pure quark EOS 
are reliable. Under this situation, we make a smooth 
interpolation of the two descriptions phenomenologically 
with a procedure simi lar to the one at high temperature 
(jAsakawa et al.l (|1997t) ): 

P = P H xf_+P Q xf + , (3) 
/ ± 4(l±tanh(^)), (4) 

where Ph and Pq are the pressure in the pure hadron 
matter and that in the pure quark matter, respectively. 
One should not confuse our interpolated pressure in 
Eq.Q with that of the hadron-quark mixed phase as- 
sociated with the first order phase transition: In our 
crossover picture, we are free from the 'phase equilibrium 
condition' such as Pq = Ph at fixed chemical potential. 
(For phenomenological attempts to interpolate H-EOS 
and Q-EOS within the conven tional picture of firs t or- 
der phase transition, see e.g. IBurgio et al.l (|2002D and 
references therein.) 

The energy density e as a function of p is obtained 
by integrating the thermodynamical relation, P = 
p 2 d{e I p) I 'dp, which guarantees the thermodynamical 
consistency for any values of p and T. To explore the 
relation between the interpolated EOS and the maxi- 
mum mass of neutron stars, we change the parameters 
(p, r) under two conditions: (i) The system is always 
thermodynamically stable dP/dp > 0, and (ii) the nor- 
mal nuclear matter is described dominantly by Ph, i.e. 
the condition p — 2Y > po. 

Shown in Fig. [5] is the pressure P as a function of 
baryon density p obtained by interpolating TNI3u H- 
EOS (the dashed line) and NJL Q-EOS (the dash-dotted 
line) with g v = Gs according to the procedure in Eq. ([J}. 
(Adopting TNI2u instead of TNI3u gives essentially the 
same final result.) As indicated by the shaded region on 
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Fig. 2. — Solid line: Pressure (P) as a function of baryon density 
p obtained by interpolating H-EOS and Q-EOS with g v = G s in the 
crossover region (the shaded region in density). Dash-dotted line: 
Q-EOS in the NJL model with g v = G S - Dashed line: H-EOS 
with the TNI3u interaction. The inset shows EOSs interpolated 
with g v /G s = 0, 1.0,1.5. 
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Fig. 3. — Solid lines: Pvs. energy-density (e) for the interpolated 
EOSs with g v /Gs = 0,1.0,1.5. Dashed line denotes hadronic 
matter with TNI3u interaction. 

the horizontal axis, crossover parameters in this figure 
are (p, T) = (3po,po), which leads to a smooth interpo- 
lation between the relatively soft H-EOS and the rela- 
tively stiff Q-EOS. Note that we have P H < P < Pq in 
the crossover region, which will never be realized in the 
standard Maxwell or Gibbs construction assuming the 
first-order transition, where P < Pjj always holds in the 
mixed phase. As will see later, the stiff EOS at (2 — 4)po 
is most important in obtaining NSs with mass greater 
than 2M Q . 

In the inset of Fig. ® we show EOSs interpolated in 
the same crossover region but with different values of 
the vector coupling (g v /Gs = 0,1.0,1.5 motivated by 
Eq.f2|)). As g v increases, EOS becomes stiffer at high 
density, while the density where the strange-quark starts 
to appear (p ~ 4po) is insensitive to g v . 

In Fig. [31 we plot the pressure P as a function of the 
energy density s for the three EOSs interpolated in the 
crossover region, (p,T) — (3po,Po)- These curves are the 
basic inputs when we solve the Tolman-Oppenheimer- 
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Fig. 4. — Solid lines: NS mass (M) vs. the central baryon 
density (p c ) obtained from interpolated EOSs between TNI3u H- 
EOS and NJL Q-EOS with g v /G s = 0,1.0,1.5. Dashed lines: 
The same quantities for H-EOS with TNI3u, AV18+TBF+AE and 
SCL3AS. The cross symbols denote the points of Af max , while 
the filled circles denote the points beyond which the strangeness 
appears. 

Volkov (TOV) equation. The case for H-EOS with 
TNI3u is also plotted by the dashed line for comparison. 
Existence of the region where P (the solid line) becomes 
larger than Pjj (the dashed line) is essential for having 
large M max as shown below. 

5. HYBRID STAR STRUCTURE 

In Fig. @J the solid lines indicate NS masses as a 
function of the central density p c obtained by solving 
the TOV equation with the EOSs given in the inset of 
Fig. [2 i.e., interpolated EOSs between TNI3u H-EOS 
and NJL Q-EOS with g v /G s = 0, 1.0, 1.5. Results with 
the hadronic EOS only (TNI3u, AV18+TBF+AS and 
SCL3AS) are also plotted for comparison by the dashed 
lines. The filled circle indicates a point beyond which 
the strangeness appears in the central core of the star. 
The cross symbols indicate the points where M max is 
achieved. 

For TNI3u H-EOS only, we have M max = 1.82M©, 

while, for interpolated EOS with g v /Gs^l, we have 
-Mmax > 2M Q . This indicates that the existence of the 
quark core inside neutron stars is compatible with the ob- 
servation of high-mass neutron star as long as (i) there is 
a smooth crossover between the hadronic phase and the 
quark phase, and (ii) there exits a strong quark correla- 
tions with repulsive nature inside the quark matter. 

In Fig. [5l the solid lines indicate the mass(M)- 
radius(i?) relation of the NSs with the interpolated EOSs 
given in the inset of Fig. [5] Results with the H-EOS 
only (TNI3u, AV18+TBF+AS and SCL3AS) are also 
plotted for comparison by the dashed lines. The gray 
band shows M = (1.97 ±0.04)M f r 1 corresponding to PSR 
J1614-2230 (pemorest et all (|2010D ). Because of the stiff 

EOS at high density for g v /Gs ~ 1, the radius R for 
0.5 M/Mq 2 takes almost a constant value, 11 — 11.5 
km tjOzel e t al. (201l])). Also, strangeness emerges only 
for massive stars with M ^ 2M . 

So far, we have taken (p,T) = (3po,Po) a s character- 
istic parameters for the crossover region. In Table I, we 
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Fig. 5. — Solid lines: M-R relation of the NSs with the in- 
terpolated EOSs between TNI3u H-EOS and NJL Q-EOS for 
g v /Gs = 0,1.0,1.5. Dashed lines: The same quantities for H- 
EOSs with TNI3u, AV18+TBF+AS and SCL3AE. The cross 
symbols denote the points where M reaches M max . The filled 
circles denote the point beyond which the strangeness appears. 
The gray band denotes M = (1.97 ±0.04) Af s for PSR J1614-2230 
UDemorest et al.l (2010).) 



TABLE 1 

M max /Mo AND pc/ P0 IN THE PARENTHESIS UNDER THE VARIATION 
OF THE CROSSOVER REGION, p AND T. 



p 


t/p 


o = 1 


T/P 


= 2 


9v=Gs 


g v =l.BG 3 


9v=Gs 


9v =1.5G s 


3po 


2.07 (5.9) 


2.18 (5.6) 






4po 


1.93 (6.7) 


2.00 (6.6) 






5po 


1.79 (7.7) 


1.83 (7.4) 


1.82 (7.4) 


1.86 (7.3) 


6po 


1.70 (8.3) 


1.70 (8.3) 


1.73 (8.0) 


1.74 (8.0) 



show how Mmax and p c depend on the choice of these pa- 
rameters for fixed values of g v . As the crossover region 
becomes lower and/or wider in baryon density, the EOS 
characteristic to the NS core (p = (2 — 4)po) becomes 
stiffer This is why M max > 2M becomes possible for 
p/po = 3 and 4. 



TABLE 2 

M max /M0 AND pc/ P0 IN THE PARENTHESIS FOR DIFFERENT 
H-EOSS WITH p = 3po AND T/po = 1. 



H-EOS 


r/po = i 


g v =o 


g v =Gs 


9 V =1-5G S 


TNI3u 


1.55 (6.8) 


2.07 (5.9) 


2.18 (5.6) 


TNI2u 


1.50 (7.4) 


2.05 (6.1) 


2.17 (5.5) 


TNI3 


1.51 (6.8) 


2.04 (6.1) 


2.16 (5.5) 


AVig+TBF+AS 


1.52 (6.8) 


2.06 (6.1) 


2.17 (5.5) 


SCL3AS 


1.55 (6.8) 


2.06 (5.9) 


2.17 (5.5) 



Now, we examine how M max depends on the different 
choice of H-EOS. Shown in Table 2 is the maximum and 
central density for different H-EOSs with p = 3po and 
T/po = 1- As long as M max is concerned, the result is 
insensitive to the choice of the hadronic equation of state 
with hyperons. 
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6. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have investigated the maximum mass 
of hybrid stars with the EOS constructed on the basis of 
the idea of hadron-quark crossover; hadronic EOS (with 
hyperons and universal three-baryon repulsion) and the 
quark EOS (with strangeness and universal two-quark 
repulsion) are interpolated to model the crossover region. 

We have found that massive hybrid stars compatible 
with M > 2Mq are possible as long as the crossover 
proceeds through relatively low density region, e.g., p = 
(3 — 4) / o u with r = po) and, at the same time, the quark 
matter is strongly interacting g v /Gs ~ 1 — 1-5. By in- 
creasing g v further, one can even obtain M max close to 
2.4 Mq. This conclusion is in remarkable contrast with 
the conventional EOS for hybrid star derived through the 
Maxwell or Gibbs construction where the resultant EOS 
becomes always softer than hadronic EOS and thereby 
leads to smaller M max . 

Although we have mainly focused on the results us- 
ing TNI3u H-EOS corresponding to k = 300MeV, very 
similar results are obtained for TNI2u as shown in Ta- 
ble 2. This means an interesting possibility that we can 
reconcile the existence of massive neutron stars (M > 
2Mq) with the ex perimental nuclear i ncompressibility 
K = (240±20)MeV (jShlomo et all (|25ol ). We note that 
TNI3, AV18+TBF+AS and SCL3AE give almost the 
same M max with that of TNI3u and TNI2u. This does 
not, however, imply the irrelevance of the universal three- 
body force: For example, the interpolated EOS with H- 
EOS (TNLSu) and Q-EOS {g v = G„) leads to strangeness 
mixing only above 4po due to the universal three-body 



force, so that the hybrid stars lighter than 1.92M do not 
have strangeness (see Fig. 4). Then, the hyperon direct 
Urea process (e.g., A — > p + e~ + D e , p + e~ — > A + i/ e ), 
which leads to extrem ely rapid neutrino-cooli ng and con- 
tradicts observations ({Takatsuk a et al.l (|2006[ )). does not 
take place except for massive NSs. 

Our massive hybrid stars with M > 2M originate 
from the transition from the soft EOS to stiff EOS in the 
density region p <~ (2 — 4)po due to the hadron-quark 
crossover. Therefore it would be important to explore 
such crossover by independent laboratory experiments 
with medium-energy heavy-ion collisions. 

Finally, we remark that the crossover region may have 
rich non-perturbative phases such as color superconduct- 
ing phases, inhomogeneous phases, t he quarkyonic phase 
and so on (jFukushima et all (|2011f )). How these struc- 
tures affect the basic conclusion of the present paper 
would be an interesting future problem to be examined. 
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